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An analytical theory of the nonlinear electromagnetic response of a two-dimensional (2D) electron 
system in the second order in the electric field amplitude is developed. The second-order polariz- 
ability and the intensity of the second harmonic signal are calculated within the self-consistent-field 
approach both for semiconductor 2D electron systems and for graphene. The second harmonic gen- 
eration in graphene is shown to be about two orders of magnitude stronger than in GaAs quantum 
wells at typical experimental parameters. Under the conditions of the 2D plasmon resonance the 
second harmonic radiation intensity is further increased by several orders of magnitude. 
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I. INTRODUCTION 

Graphene is a recently discovered^ purely two-dimensional (2D) material consisting of a monolayer of sp 2 -bondcd 
carbon atoms arranged in a hexagonal lattice. Electrons and holes in graphene are massless Dirac fermions and this 
leads to a variety of interesting and unusual electrical and optical properties of this material—. It promises a lot of 
applications in electronics, optics and other areas^- 7 - . 

It has been predicted^ that the unusual linear energy dispersion of the charge carriers should lead to a strongly non- 
linear electromagnetic response of graphene: the irradiation by electromagnetic waves should stimulate the emission 
of higher frequency harmonics from graphene. The theory of the nonlinear electromagnetic response of graphene has 
been further developed in Refsi^r— . Experimentally the higher harmonics generation and frequency mixing effects 
have been observed in Refst^ 7 - - — . 

In theoretical papers^— only the normal incidence of radiation on a uniform graphene layer has been studied. The 
experimental demonstration 2 ™ of the third-order emission of radiation at the frequency 2lui — UJ2 at the bichromatic 
irradiation by the frequencies w\ and ui2 has confirmed that graphene manifests the nonlinear properties and that its 
third-order effective nonlinear susceptibility is much higher than in a number of other materia ra 15 i 16 i 20 . However, the 
intensity of the emitted signal /2wi±u 2 ^ ^2i-^2 is proportional to the third power of the intensities of the incident 
wavesSS therefore to observe the third-order nonlinear effects one needs quite powerful sources of radiation. 

Substantially stronger nonlinear effects could be expected in the second order in the radiation electric field. The 
second-order effects, e.g. the second harmonic generation, are proportional to the second power of the incident wave 
intensities, oc / 2 . However graphene is a centrosymmetric material, therefore at the normal incidence of radiation 
the second-order effects are forbidden by symmetry. 

The symmetry arguments do not hinder the observation of the second-order effects at the oblique incidence of 
radiation on the 2D electron layer. If the incident wave has the wavevector component q parallel to the plane of 
the 2D layer, one could observe much stronger second-harmonic radiation as compared to the third order effects 20 . 
Moreover, at the oblique incidence of radiation one can resonantly excite the 2D plasma waves in the system (e.g. in 
the attenuated total reflection geometry) which would lead to the resonant enhancement of the higher harmonics 2 ^. 

In this paper we theoretically study the second-order nonlinear electromagnetic response of two-dimensional electron 
systems, including both graphene and conventional 2D structures with the parabolic electron energy dispersion. We 
calculate the second-order polarizability of 2D electrons (Section [H]) and show that in graphene it is at least one 
order of magnitude larger than in typical semiconductor structures (e.g. in GaAs/AlGaAs quantum wells). Then 
we calculate the self-consistent response of the system to an external radiation taking into account the 2D plasmon 
excitation (Section IIII[) and show that the intensity of the second harmonic signal can be further increased by several 
orders of magnitude. In Section ITVl we summarize our results. 
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II. RESPONSE EQUATIONS 
A. General solution of the quantum kinetic equation 

We consider a 2D electron system with the Hamiltonian Hq, where 

H \\) = E X \\), (1) 

A is a set of quantum numbers, E\ and |A) are the eigenenergies and eigenfunctions of the system. Under the action 
of the electric field described by the scalar potential 

0(r,t) = ^ e ^ r ^ t )+c.c. (2) 

the state of the system is perturbed (here c.c. means the complex conjugate). Its response to the potential (f5]) is 
determined by the Liouvillc equation 

i?i^ = [H,p] = [H + H 1 ,p] (3) 

where p is the density matrix and H± = —e<f>{r,t). Our goal is to calculate the charge density fluctuations 

p(r, t) = p^S**-*) + P2^e 2 ^ r ~^ + c.c. (4) 



in the first and second orders in the potential amplitudes (/> qtlJ and the corresponding first- and second-order polariz- 

Pquj — '-'qJiqu'/'qu I (5) 



abilities 42 ;qw and a^q^qi^qw defined as 



(6) 

In the absence of the perturbation H± the density matrix po satisfies the equation 

Po|A> = / A |A), (7) 

where fx — f(E\) is the Fermi distribution function. Expanding p in powers of the electric potential, p = po + p\ + 
f>2 + • ■ ., and calculating the charge density fluctuations — eSp [S(r — Tq) {pi + p2 + ■ ■ ■)] we get 

A A' 



a 



sl— = 4g £ ,-tTZU,o i:' A ' e,q, ' A ">< v, ' e ' qr ' y > 

/ A' — f\" f\" — fx \ 



E X ' - Ex" + Tiuj + iO E x >> - E x + hco + iO 



(9) 



where S is the sample area. The first-order polarizability (JSj> is proportional to the polarization operator II(q, lu) 

(see22,), aqj : qcj = — e 2 II(q, uj). For the conventional 2D electron gas (with the parabolic electron energy dispersion) 
the linear polarizability ([S]) has been calculated in Ref.— , for 2D electrons in graphene (with the linear energy 
dispersion) it has been done in Refs . 23 i 24 . 

We will apply the general formulas ©, © to the conventional 2D electron systems in semiconductor heterostruc- 
tures and to graphene. In the former case the spectrum of 2D electrons is parabolic, in the latter case it is linear. In 
both cases the set of quantum numbers |A) = |Zkcr) consists of the subband index I, the wavevector k and the spin 
a. To specify the general expressions ([5]), © we consider the long- wavelength limit, which is quantitatively described 
by the conditions 



q max{fci?, fcy}, q *C w/ max{w^ , vt} 



(10) 
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in semiconductor 2D electron systems and the conditions 

gCmaxfip,^}, q^iui/vF (11) 

in graphene. Here kp and vf are the Fermi wave- vector and Fermi velocity respectively kp and vp are the thermal 
wave-vector and velocity respectively. In the gas with the parabolic dispersion kp = ^/2mEp /ti, kp — \/2mT/h, 
vp.T = frkp^p/m, where m is the effective electron mass, T is the temperature and the Fermi energy Ep is counted 
from the bottom of the parabolic band. In the 2D gas with the linear energy dispersion (in graphene) kp = y/\JI\/hvF, 
kp = Vr/hvp, where fi is the chemical potential counted from the Dirac point (/z can be positive and negative) and 
Vf is the Fermi velocity. At typical experimental parameters the conditions (| 10[) and (jlip restrict the wavevector q 
by the values ~ 10 6 cm -1 which is sufficient for the description of most experiments. 

Under the conditions (jTUJ) - (fTTj) the general expressions © and can be substantially simplified and we get the 
first and second order polarizabilities of the 2D electron gas in the form 

(i) ^ e 2 g a g/9 9b \^ ( df ik \ dE lu 



(2) _ 3e 3 'q a q/3qjqs g s \ - d/ik gg[k d 2 E i ]i 

*2qS 

where g s = 2 is the spin degeneracy. 



., , k dk a dk p dk^dks' [ ' 



B. 2D electron gas with the parabolic energy dispersion 



Let us apply now the obtained formulas (|T2l and (flUl) to the conventional 2D electron gas with the parabolic energy 
dispersion. In this case there is only one energy subband (/ = 1), 

E lk = E k = ^ (14) 
Am 

and the wave functions are plane waves. Substituting (fT4|) in Eqs. (fl~2j) and (fT3|) we get 

(1) _ n s e 2 g 2 

TOW 2 ^ 15 ) 

and 

(2 ) _ 3n s e 3 q 4 

a 2q2u;;q^qw ~ 2TO 2 CJ 4 ' 

In semiconductor 2D electron systems both o^ 1 ) and are proportional to the 2D electron gas density n s . The 
result (1151) has been obtained in Refj22. 



C. 2D electron gas with the linear energy dispersion (graphene) 



In graphene the spectrum of electrons can be found in the tight-binding approximation^. It consists of two energy 
bands, the wave functions are described by Bloch functions and the energy of the electrons is 

E^ = (-l)'t|5k|, 1=1,2, (17) 

where t is the transfer integral, 5k is a complex function defined as 

5 k = 1 + e lkai + e ika2 = 1 + 2 cos(k x a/2)e zV3k y a/2 , (18) 

ai = o(l/2, n/3/2), a 2 = o(-l/2, V5/2) are the basis vectors of the graphene hexagonal lattice and a = |ai| = |a.2 1 is 
the lattice constant. In graphene t ps 3 eV and a = 2.46 A. 



4 



The formulas ([H])-© and (fl"2"|) - (fT3]) are valid for the full graphene energy dispersion, i.e. the integration in these 
formulas is performed over the whole Brillouin zone. Under the real experimental conditions, when the chemical 
potential /i satisfies the condition \/j,\ <C t, the main contribution to the integrals (fl~2|) and (fl~3f is given by the vicinity 
of two Dirac points K x = (27r/a)(l/3, 1/V3), K2 = (27r/a)(2/3, 0), where the function 5k vanishes, <5>Ki = <Sk 2 = 
and the energy ()17|) can be approximated by linear functions 



E lk = (-l)W|k| = (-l) l hv F \k - K v \, v = 1,2, 



(19) 



where v is the valley index and k = k — K„ is the electron wavevector counted from the corresponding Dirac points. 
The velocity parameter vp here (the Fermi velocity) is related to the transfer integral and the lattice constant, 
vp = in graphene vf ~ 10 s cm/s. Omitting below the tilde over the wavevector k and calculating the 

derivatives 



BE, 



Zk 



dk a ^ ^ F k ' dk a dkfj 
we get the first order polarizability in the form 



{-i) l hv F - 



1 a p 



fc 3 



-. 

where g v — 2 is the valley degeneracy and 

F(a) = 



4irh 2 u 2 



F 



2T 



F 



JAM 
2T 



xdx 



o cosh (x — a) 



2a, a > 1 
ln2, a = 



(20) 



(21) 



(22) 



If the temperature is low as compared to the chemical potential, T <C we get from (|2ip the result obtained in 
Refs3^, 



e 2 g s g v q 2 \fi\ 
4tt^ 2 w 2 : 



T « 



In the opposite case <C T one has 



e 2 gs9vq 2 T 

ATTh 2 UJ 2 



21n2 



(—)' 

\2TJ 



ImI « r. 



Now consider the second order polarizability of graphene. Using Eqs. (|TH)) and (|2U)) we get 



(2) 

2q2cj;qcj,qcj 



(2) 

The second-order graphene polarizability a^q^ qu qw ^ s an °dd function of the chemical potential and does not 
depend on the electron/hole density at > T. This is the direct consequence of the linear energy dispersion (|19p 



3e 3 g s g v q 4 v F ^ jj_ 



(23) 



(24) 



(25) 



and essentially differs from the case of the conventional 2D electron gas (jTHJ) for which cx n s . The q- and 
w-dependencies of in the linear- and parabolic-spectrum cases are the same, cx q 4 /uj 4 . 

Let us compare the second-order polarizability of graphene (|2"5|) with that of a conventional 2D electron gas (|16p , 
e.g. for typical parameters of 2D electrons in a GaAs/AlGaAs quantum well. Assuming that m = 0.067mo and 
n s = 3 x 10 11 cm~ 2 for GaAs (mo is the mass of free electron) and \fi\ ^> T for graphene we get 



* (2) „ 

graphene 

(2) 
a GaAs 



10. 



(26) 



III. SELF-CONSISTENT FIELD 



Equations §5§ and ([6]) determine the first- and second order response of the 2D electron gas to the electric field 
really acting on the electrons. Consider now the experimentally relevant formulation of the problem when the system 
responds to the external field 4>^J ■ Using the self-consistent field concept we solve, first, the linear response and then 
the second-order response problem. 
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A. Linear response 



Consider a 2D electron system under the action of an external electric potential (j) ext (r , z , t) — (p^ x J(z)e' lcl ' r ^' lut . 
In the first-order in the external field amplitude the resulting 2D charge density will also contain the qw-harmonic 
p(r,t) — /0 qtlJ e 4q ' r_lw *. The density fluctuation creates, in its turn, the induced potential 4>™ d (r, z,t) = (j) l ^e lcl ' r ~ %ut 
determined by Poisson equation 

A(j) md (r, z, t) = -47rp(r, t)S(z) (27) 

and given by 



< d (* = 0) = ^ d = ^ q ,.. (2Si 



The density pq U here is determined by the response equation ([5]) in which the really acting on the electrons potential 
should stay in the right-hand side. This is not the external field but the total electric field <^fj = ^qw + ^qu produced 
both by the external charges and by the 2D electrons themselves. Then we get 



bqu — fiqu ^qu] ~ Pqu — a q!;quj ( l ) qu (29) 



and the known linear-response formula 



with the dielectric function 



<b ext 

*£ = 7&TT. (3°) 



e(q,w)' 



e(q,w) = l-^a« qlir (31) 

If the wavevector q and the frequency uj of the external wave satisfy the condition 

e(q,w)=0, (32) 

one gets a resonance in Eq. (|30|) . This resonance corresponds to the excitation of the eigen collective modes of the 
system - the 2D plasmons. 

Consider the 2D electron gas with the parabolic energy dispersion. Substituting the linear polarizability (|15j) into 
Eq. (f3"Tj) we get from (f3"2"|) the known spectrum of the 2D plasmon 

27T7J 2 

u* = u$(q) = —?—q (33) 

first obtained in Ref^ (we ignore the dielectric constant of the surrounding medium; if the 2D gas is immersed in 
the insulator with the dielectric constant k, e 2 here should be replaced by e 2 /n). 

In the case of graphene, the spectrum of the 2D plasmons follows from Eqs. (|2"T|) . (|31l) and (1321) : 



2 2 1 \ e 9s9vT 

bj = 0J n (q) — t. — 

pK ' 2Ti 2 



' (34) 



F ,4 ) +F _M 



2T \ 2T 



In the limit \fi\ 3> T this gives the result 



%, s e 2 g s g v \fi\ 

Up(q) = — ^2 — ? ( 35 ) 



obtained in Refs . 23 ' 24 . In the opposite case \fi\ <cTwe get 

2/ \ e 2 g s g v T \n2 

Wpfe) = p <?■ (36) 

The 2D plasmon problem in the regime /i = has been considered in Ref.— . The result reported in 2 ^. differs from the 
correct formula (|3"6")l by a factor of Air. 
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B. Second order self-consistent response 

Let us now consider the second order response to the external potential <j)^e lcl ' T ~' lu!t + c.c. The induced and total 
potential will now contain the frequency harmonics ±(qw) and ±2(qw). The self-consistent charge density then reads 

n(r t\ — /?W Atot iq-r-iut , (1) -tot i2q-r-i2wi , (2) rtotitoti2q-r-i2wt, (o 7 \ 

where the first two terms correspond to the linear response to the first and second harmonics and the third term - 
to the second-order response to the first (qw) harmonic of the total potential. The complex conjugate terms describe 
the negative (qw) harmonics. The second order response to 0q^<^*f( qw ) vanishes. 

The Fourier harmonics of the induced potential follow from Eq. ([37| and Poisson equation (|27[) : 

6 Jr t\ - — tyi 1 ) jJotiq-T-iut , ^ (1) -tot i2qr-i2wt , ^ (2) -tot -tot i2q-r- l 2 W t , 

nn<il r !^ — ^ "quiquV'qu 6 T £^ a 2q2w;2q2tjV , 2q2w e ' 2g "2q2t l ;;qw,qwS y qwS t 'qLj e; ' C.C 

lind iq-r—iuit , lind z2q-r— z2ut , -oo\ 

= 0qc;e q +<P 2 q2y e + CC. (38) 

Equating now the amplitudes of the first-order harmonic (qw) we get from here Eqs. (|29|) and pO|) . Equating the 
coefficients at the (2q, 2w)-harmonic and taking into account that the second harmonic component is absent in the 
external potential, 02q2w = 0, we get 

lind itot _ ^ A 1 ) Jfot , n (2) itot itot (oq\ 

S y 2q2w — V 9 2q2u _ r, a 2q2w;2q2wS y 2q2w T a 2q2w;qw,qw ( ?V l jY ; qw l JJ 7 

q 

and finally 

_,(2) (2) 

-tot " 2q2w;qw,qw , tot , tot " "2q2w;qw,qw text text { Ari \ 

= g C (2q,2 W ) = ^ e (2q,2.)[ £ (q,.)] 2 ^ ^ • (40) 

The formula (|4T))) . together with ([T5| and (|23|). represents the main result of this work. One sees that the amplitude of 
the second-harmonic potential is resonantly enhanced at the frequency ui — uj p (q) [the second-order pole corresponding 
to the zero of e(q, u)] and at the frequency u> — ui p (q)/y/2 [the first-order pole corresponding to the zero of e(2q, 2ui)]. 
These resonances allows one to get a huge enhancement of the second harmonic radiation intensity. 



C. Estimates of the second harmonic radiation intensity 



Let us estimate the intensity of the second harmonic signal. Assume that the external potential is 

(j} ext (r, t) = cf> cos(q • r — wt), 
so that q ^' = (j>o/2. Then the total potential at the frequency 2ui obtained from Eq. (|40)) reads 



6 tot , t) _ Wo (2) 

r2q2cj V ' / — 2g a 2q2w;qoj,qw 



(2(q-r-w*)) 



( W 2_^y 



where we have introduced the momentum scattering rate 7 in the dielectric function 

<4(q) 



e(q,w) = 1 



(41) 



(42) 



(43) 



Text 
qcj 



to remove unphysical divergencies of the plasma resonances. Introducing now the intensity of the incident J q 
cEq/8tt and of the second-harmonic wave /| q 2^ ~ cEf ot /8ir, where Eq and E tot are the fields corresponding to the 
potentials (|41l) and (1421) respectively, we get the following results: 



1. In the 2D electron gas with the parabolic electron energy dispersion (conventional semiconductor structures) 



18ir 3 n 2 s e 6 q 4 



2 ext. • 



(44) 
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2. In the 2D electron gas with the linear energy dispersion (graphene) 



I, 



graphene 
tot 



=^tanh^(A) 



8h 4 c 



U</> 



(w 2 -W2(g)) 2 + W 2 7 2 



-r 

2 ext " 



(45) 



The ratio of the intensities (|4"5")) and (l4"4l is proportional to the squared ratio of the polarizabilities a*- 2 - 1 . For the 
same parameters that have been used in Eq. (|26|) one gets 



jgraphen 
1 tot 

jGaAs 

1 tot 



graphene 

~ (2) 

a GaAs , 



100. 



(46) 



The frequency dependence is the same in both graphene and semiconductor cases and is shown in Figure [TJ The 
intensity versus frequency curve has a huge resonance at the frequency u> ~ uj p (q) and a weaker one at the frequency 
u) ~ u) p (q)/\/2. The chemical potential and temperature dependence of jaraphene ^ s \ lown m Figure [2] 



If uj = oj p (q) (the main resonance maximum) the ratio If™ phene / I e xt can be presented as (at S> T) 



jgraphene 
L tot 



9TTe 6 v 4 F q 4 Ie X t x ^p(g) 
2h A cwl{q) X 7 4 



3^o 



u p (q) 



(47) 



In this formula Eq is the electric field of the external incident electromagnetic wave and ekp is the internal electric 
field in the 2D system (the field produced by an electron at the average inter-electron distance kp 1 ). The ratio Eo/ekp 
in the first brackets is therefore typically very small, Eo/ekp <C 1. The second factor u p (q)/^i is the quality factor of 
the 2D plasmon resonance, which can be very large in the high-quality samples. This may, at least partly, compensate 
the smallness of the first factor and substantially facilitate the observation of the second harmonic generation. 



IV. SUMMARY 



We have presented the self-consistent analytical theory of the second harmonic generation in two-dimensional 
electron systems. The theory is applicable to semiconductor structures with the parabolic and graphene with the 
linear electron energy dispersion. We have shown that the intensity of the second harmonic is about two orders of 
magnitude larger in graphene than in typical semiconductor structures. Under the conditions of the 2D plasmon 
resonance the intensity of the second harmonic can be enhanced by several orders of magnitude. 

The frequency of the 2D plasmons in graphene lies in the terahertz ranged—. The discussed phenomena can be 
used for creation of novel devices (frequency multipliers, mixers, lasers) operating in this technologically important 
part of the electromagnetic wave spectrum. 
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FIG. 1: The intensity of the second harmonic radiation in semiconductor structures and in graphene as a function of the 
frequency u)/u p (q). The parameter of the curves is r y/ui p {q). 




FIG. 2: The intensity of the second harmonic radiation in graphene as a function of the chemical potential and temperature 
M/2T. 
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